Diagonal patterns and chevron effect in intersecting traffic flows 
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We study a lattice model of two perpendicular intersecting flows of pedestrians represented by 
hard core particles of two types, eastbound ('£') and northbound ('A/ - '). Each flow takes place 
on a strip of width M so that the intersection is an M x M square lattice. In experiment and 
simulation there occurs on this square spontaneous formation of a diagonal pattern of alternating £ 
and N particles. By a linear stability analysis of the corresponding mean-field equations we point 
out the origin of this pattern formation phenomenon. A refined investigation reveals that for large 
enough M the pattern actually consists of chevrons rather than straight diagonals. We provide an 
explanation for this chevron effect. 
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Pedestrian motion in dense environments is of 
much theoretical and practical interest [H, Q ■ 
Instances of applications are shopping streets, 
waiting lines, crowds that enter or leave a room, con- 
verge to a stadium, participate in a demonstration, and 
so on. Simplified models help us understand the behav- 
ior of individuals under such circumstances as well as the 
collective behavior that results from it. They may also 
have more general applicability than the question that 
initially prompted their formulation. 

The present investigation was motivated by the ex- 
perimental observation 0] of an instability that occurs 
at the intersection of two perpendicular unidirectional 
flows of pedestrians: Walkers segregate into a pattern 
of approximate diagonals of walkers of the same type. 
In a variety of models simulating crossing flows (agent 
based HQ, PDEs 13. BML 0) it; has been remarked 
that such patterns occur, but the formation mechanism 
has not been systematically studied. In this paper we 
tackle this question with combined numerical and theo- 
retical approaches. First, we explain the pattern forma- 
tion instability both for closed and open systems. More- 
over, we show that superposed on the instability there is 
a subtle 'chevron effect', whose origin we explain. This 
effect, which so far has not been observed, may be ex- 
pected to be visible under favorable circumstances. 

We model a street of width M as a set of M parallel 
one-dimensional lattices or 'lanes' Two intersecting 
one-way streets lead to the geometry of Fig.[TJ which has 
the (1,1) diagonal through the origin as an axis of sym- 
metry. Eastbound (or '£') and northbound (or 'A/') par- 
ticles are injected into each lane at some large distance L 
from the intersection square with a probability a per time 
step. This probability determines the incoming current 
J(cv) in each lane For moving the particles we choose 



for convenience the frozen shuffle update It will ap- 
pear that the phenomena of interest are independent of 
the details of the update mechanism. Under frozen shuf- 
fle update, during a unit interval all particles are visited 
in a fixed sequence and each one executes a move un- 
less its target site is occupied; that is, the particles move 
ballistically unless blocked. Particles entering the system 
are randomly inserted in the update sequence and those 
exiting are deleted from it. Blocking may happen to a 
particle either in the intersection square itself or in the 
street segments leading up to it. In the simulations we 
chose L (see Fig.QJ larger than the lengths of any waiting 
lines observed at the entrance, so that effectively L = oo. 
We therefore have a simple model depending only on the 
two parameters a and M, whose interest as an example 
of a driven noncquilibrium system goes, we believe, be- 
yond the present application. Simulations were carried 
out for intersecting streets of widths up to M = 640. The 
pattern formation studied here appears in the so-called 
free flow phase (a < 0.10 for the M values concerned), 
in which in each lane the current exiting the intersection 
is equal to J (a). We will stay away from higher values of 
a, where jamming transitions are known to occur [jlfioj]. 

For t > particles start entering the initially empty in- 
tersection square. After a number of time steps typically 
no larger than a few times the linear lattice size M, the 
occupation of the intersection square reaches a station- 
ary state. We observed that this state has the following 
properties. 

(i) There is an a dependent penetration length £(a), 
such that for M < £ the occupation of the intersection 
square appears disordered to the eye. For a = 0.09 the 
penetration length equals £(a) « 50; for a — > it di- 
verges. 
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FIG. 1: Intersection of two one-way streets of width M. The 
blue particles (►) move eastward and the orange particles 
(a) northward. The parameter a determines the particle in- 
jection rate. The region bordered by the heavy solid line is 
the 'intersection square'. 




FIG. 2: Typical configuration of the intersection square in the 
stationary state for M = 640 and a — 0.09. The orange par- 
ticles arrive from the left and the blue ones from the bottom. 
Between the lower-left and the upper-right, the particles self- 
organize to form a diagonal pattern. The dashed white lines 
delimit different regions discussed in the text. 



(ii) For M > £ alternating stripes of £ and Af par- 
ticles parallel to the (1,-1) direction begin to be dis- 
tinguishable. For Al ^> £ they are clearly visible, as 
may be seen from Fig. [2 where M = 640 and a = 0.09. 
Although the striped pattern never becomes fully reg- 
ular, its wavelength is typically in the range from 5 to 
15 lattice distances. This organization into stripes de- 
creases the probability for a particle to be blocked below 
its value for a random particle distribution, and therefore 
increases the particles' average velocity. 

(hi) Closer examination of Fig. [5] reveals an effect just 
barely visible to the eye: the angle 6 of the striped pat- 
tern (as measured clockwise from the west) is not exactly 
equal to 45° but deviates from it by an amount A8(r), 
where r = (x, y) denotes the position in space. Although 
small (of the order of a degree), the deviation A6*(r) may 
be measured (l2j unambiguously in the neighborhood of 
any site r. It is negative above the axis of symmetry and 
positive below it. 

(iv) In the upper (lower) triangular region delimited 
by the dashed white lines in Fig. [21 the deviation A9(r) 
is close to a constant =fA#o- This confers to the stripes 
the appearance of chevrons, and we will speak of the 
'chevron effect'. The tips of the chevrons, located in the 
transition zone between the two white lines, are rounded. 
Empirically we find A8o(a) ~ ca with c « 15°. 

We will now first explain the stripe formation instabil- 



ity and then return to the chevron effect. 

Let us denote by p\ "^(r) the particle densities of the 
£,J\f particles on site r after t time steps. In an exact de- 
scription of the time evolution the p\ ' (r) would couple 
higher order correlations. Breaking these, we postulate 
the mean-field equations 

pf +1 (r) = [1 - p?(tM (r " e«) + (r + e x )pf (r), 
p? +1 {v) = [1 - p £ t {v)]pf{r e y ) + pf (r + e y )p? (r), 

(1) 

where e x y are basis vectors. Eqs. ([T]) take into account 
the interactions between particles of different type (E and 
J\f) but neglect those between particles of the same type. 
This may be justified by the fact [l2[ that in the low 
density limit the frequency of £/N and Af/£ blocking 
events is ~ ~p 2 , whereas that of same-particle blocking is 
~ p 3 , where p is the average density of each of the two 
particle types in the intersection. 

As a first step we impose toroidal boundary conditions 
on ([1]) and perform a linear stability analysis around a 
uniform initial state. Let pf"^(r) = ~p+5pf ,A ^(r), where p 
now replaces a as the control parameter and the 5pf^(r) 
arc small. In terms of the Fourier transforms Spf ,JV "(q) 
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with q = (q x ,q y ) the linearized equations read 
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1 e^R qy ) 



6p?(q) 



(2) 

where R q = 1 + p(e 19 — 1). The larger eigenvalue of the 
2x2 matrix in (J5|) has an absolute value whose maxi- 
mum, located on the diagonal q x = q y , exceeds unity; it 
is associated with an unstable mode traveling in the (1,1) 
direction and having the wavelength A max = 2-k /q max = 
6[1 — (\/3/7r)p]- This calculation therefore explains the 
formation of a diagonal striped pattern as a consequence 
of unstable Fourier modes present in the initial state. 

The nonlinear regime is outside of the reach of this 
analysis; however, numerical solution of the nonlinear 
equations ([T]), still with toroidal boundary conditions, 
shows that the solution tends to a stationary state con- 
sisting of stripes with alternatingly only pf (r) ^ or only 
Pj^(r) 7^ 0, and separated by unoccupied sites in such a 
way that the nonlinear terms in ([T]) vanish at all times. 
The density profiles p £ ,A ^(r) then advance unaltered east- 
ward and northward around the torus, respectively, with 
unit velocity. 

In the original crossing street system the interaction 
square has open west and south boundaries that are sub- 
ject to time-dependent boundary noise. It is therefore 
necessary to show that this noise excites unstable modes 
analogous to those found above for the torus. For the 
open system we consider the nonlinear equations ([T]) with 
the stochastic boundary conditions 



p £ (0,k)= V £ (k), 
pf(M + l,fc) =0, 



p?(k,0) = rj?(k), 



Mi 



pf (k,M + l) = 0, 



(3) 



for all k = 1, 2, ... , M, where the j]\' N (k) are i.i.d. ran- 
dom variables of average rj and which we take from a 
uniform distribution on [iyy, |?y]. Here rj replaces a as 
the control parameter. The initial condition is arbitrary. 

After linearizing these equations we found, at the cost 
of considerable algebra, their Green function 13| and es- 
tablished the existence of unstable modes propagating in 
the (1,1) direction. This completes our explanation of 
the origin of the pattern formation instability. 

We now turn to the chevron effect. In the linear sta- 
bility analysis of (JT]) with the stochastic boundary condi- 
tions ([3]) this effect does not appear. However, numerical 
resolution of the full nonlinear equations ([T]) with ([3]) 
leads to the same chevron effect as simulation of the par- 
ticle system. This is one more strong indication that it 
must exist for a wide class of time evolution rules, as long 
as these incorporate the basic ingredient of repulsion be- 
tween two unidircctionally and perpendicularly traveling 



particle types. 

We now present further observations made on our sim- 
ulation of the particle system, which will help us explain 
the chevron effect. 

(v) The stripes of the £ and those of the Af particles 
move almost without any mutual penetration. 

This is compatible with an inclination 9 of the stripe 
direction only if tan#(r) ~ (r) / v £ (r) , where v £ '^ is 
the time-averaged velocity of the £, Af particles on site r. 
Now the boundary conditions impose the same station- 
ary current J in each lane, so that J = p £ (r)v £ (r) = 
p^ (t)v^ (r) with p £ -^(r) the time-averaged densities. 
Combined with the preceding equation for tan#(r) this 
yields 



tan0(r) = p £ {r)/ P H {r) 



(4) 



It follows from (|4]) that A6*(r) ^ requires p £ {v) and 
p^{r) to be different. We have verified that relation (|4| 
is satisfied in good approximation both in the simula- 
tions of the particle system and in the numerical solution 
of (JTJ) with J3]), in the two triangular regions as well as 
in the transition zone between them. We remark paren- 
thetically that since A8 is small, the difference between 
the two densities is small as well. Indeed, a particle 
simulation for M = 500 and with current J(a) = 0.06 
yields in the upper triangular region p £ w 0.0606 and 
k, 0.0621. This establishes the general consistency of 
the picture, still without explaining it. 

(vi) Next we refer to the zoom, shown in Fig.[3j on an 
area located in the upper triangular region of Fig. [2] We 
observe that in this region the stripes of the Af particles 
are dense, whereas those of the £ particles are sparse. A 
consequence, visible even if barely so, is that the upper 
triangular region in Fig. [2] looks bluish and the lower tri- 
angle more orange-like. The asymmetry observed here 
between the two particle types will offer the clue to un- 
derstand the chevron effect. 

The core of the problem is to show the existence of 
modes of propagation having stripe angles different from 
45°. Let us consider what happens near the entrance 
boundary of the £ particles, that is, for x w £ but y ^> £. 
Near this boundary the £ particles (►), having entered 
the intersection square randomly, fill the space offered to 
them between the Af stripes ( ) also largely randomly. 
This suggests to consider the special class of Af stripes of 
which an example is shown in Fig-Ht 1 - The stripes consist 
of straight segments at an angle of 45°, connected by 
'kinks' such as the one that occurs in Fig.[4ji at the level 
of particle B, and that is associated with the presence 
of the eastbound particle A. The other £ particles in 
Fig. 2^. occupy random positions. In order to facilitate 
the analysis we now invoke the irrelevance of the details 
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of the time evolution rules. This allows us to employ in 
the present argument the alternated parallel update: at 
each time step first all £ particles move simultaneously 
(unless blocked) and then all unblocked M particles do 
so. 

When this update procedure is applied to the config- 
uration of Fig.[4ji, we see that during the (t + l)th time 
step none of the eastbound particles is blocked; in par- 
ticular, A and C move and block B and D, respectively 
As a result the kink associated with A moves one step 
to the right along the northbound stripe, and C creates 
a new kink at the beginning of the stripe, at the level of 
particle D. (see Fig. |4Jd). In the next time steps parti- 
cles A and C will both travel from left to right along the 
stripe, each of them taking its associated kink along, and 
the connected structure of the stripe will be preserved. 

If the set of kinks has a linear density pkink along 
the stripe, the average stripe angle 8 will be given by 
tan 9 = 1 — Pkink- Since pkink also represents the fraction 
of blocked moves of the N particles, the stripe's speed 
will be lA = 1 — pkink- Hence we have demonstrated 
the most distinctive ingredient of the chevron effect: the 
existence of a nonlinear mode consisting of a stripe with 
an average slope different from 45°, that propagates at 
an average speed < 1. 

We must expect similar modes to be present for a wide 
class of time evolution rules, including those of the origi- 
nal particle system with frozen shuffle update, and those 
of Eqs. |T]); for both of these the explicit analysis would, 
however, be much more difficult. Of course quantitatives 
features such as the exact value of AO or the wavelength 
of the pattern will depend on the specificity of the time 
evolution rules. 

Returning to the example of Fig.0] we note that a 
uniformly random spatial distribution of the £ particles 
would lead to pkink = p £ , and that v £ = 1. Combining 
these approximations we are led to a fully explicit ex- 
pression for the angle, namely tan 8 = 1 — </, in which 
J (a) is known 0. Since a correlated distribution of the 
£ particles would lead to a lower pkink, we expect that 
this formula, while giving the correct order of magni- 
tude, overestimates the inclination. Indeed, for J — 0.06 
it yields A&o = 1.8°, to be compared to 0.7° obtained 
via Eq. Q and 0.9° from direct measurement, both in a 
frozen shuffle update simulation. 

We have presented a set of simulation results, numer- 
ical work, and analytic work that coherently bring out 
a pattern formation instability accompanied by a subtle 
but unmistakable 'chevron' effect, quantified by an an- 
gular deviation A0(r) of the pattern. The effect may 
be measured in the simulation and also calculated nu- 
merically. Determining the precise values of the time- 




FIG. 3: Zoom on a region of Fig. [2] To emphasize their dif- 
ference the northbound particles ( ) have been indicated by 
larger symbols than the eastbound particles (►). The black 
solid lines are at an angle of 45°. The nonzero distances 
Ai, A 2 , and A 3 show that locally A0(r) < 0. 
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FIG. 4: Mechanism causing the deviation A9 of a stripe of 
northbound particles ( ) in the upper triangular region. 



averaged density profiles p £,AA (r) of the two intersecting 
particle flows, and thence of the preferred angle of incli- 
nation A0(r), is still an open question. Its solution might 
require a hydrodynamic approach to this problem, which 
however is beyond the scope of this communication. In 
any case, we expect the chevron effect to be generic for 
a wide class of transport problems, including those of 
Refs. 0, 0, 0]> that share the features of intersecting uni- 
directional flows. Because of its smallness, observing the 
effect in an experiment or in a real-life traffic situation 
will not be easy, but must be anticipated. 
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